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Abstract 

The knowledge of capacitance in semiconductor micro-strip detectors is impor- 
tant for a correct design, simulation and understanding of the detectors. 
Analytical approaches can efhciently complement numerical methods providing 
quick results in the design phase. 

The conformal mapping method has proved to be the most effective analytical 
approach providing many realistic models 

In this paper improved analytical results are presented and compared with ex- 
perimental data. 

The excellent agreement between predictions and measurements confirms the 
relevance of this approach to modeling realistic detectors. 
Key words: Elliptic functions, Conformal mappings. Micro-strip detectors 
PACS: 33E05, 30C20 

1. The role of capacitances in micro-strip detectors 

A micro-strip detector can be modeled as a circuit: the energy release is 
a current generator and the resulting signal propagates in a network of ca- 
pacitances and resistances. The nodes of the network are the strips and the 
backplane. A realistic model might incorporate also the read out electronics 
The capacitive elements of the circuital model are shown in Fig[T] Cg is the strip 
to ground capacitance, while C„ is the interstrip capacitance n strips away. 
The relation between the output signals and the energy release is mediated by 
the detector capacitances. These are difficult to calculate because of the com- 
plex multi-electrode geometry. 

A numerical solution of the Poisson equation within the detector is possible 
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Figure 1: Capacitive circuital model of a micro-strip detector section. Every 4th strip is read 
out, employing the principle of capacitive charge division. 



but is of limited use for the design because it has to be repeated for every 
configuration. Therefore it is interesting to investigate approximate analytical 
solutions. 



2. Capacitance calculations v^rith conformal mapping 

2.1. Conformal mapping 

A powerful approach to calculate electrostatic fields and capacitances in com- 
plex geometry with 2D symmetry is the use of conformal mappings 
They are angle preserving 1 to 1 functions between regions of the complex plane. 
This property is equivalent to the requirement that they are differential as com- 
plex functions 

An electrostatic field satisfying the Laplace equation in a region with 2D symme- 
try can be described by a complex potential function F{z — x -\- iy) = t{x, y) — 
u{x,y) + iv{x,y), where u{x,y) is the potential function and v{x,y) the flux 
function defined such that the curves with constant v{x,y) and with constant 



2 



u{x,y) are orthogonal. Both v{x,y) and u{x,y) satisfy the Laplace equation 
and F is differential and conformal. 

Any conformal function preserves potential and flux functions and the capaci- 
tances between corresponding conductors. 

This property is the key for using conformal mapping in capacitance calculation. 
A region with a non-trivial distribution of conductors (and/or boundary con- 
ditions) is mapped onto a region (possibly through intermediate configurations 
and exploiting symmetry) where the capacitance is known, e.g. parallel plates. 

2.2. Schwartz- Christoff el transformations 

The class of conformal mappings relevant for our application are the Schwartz- 
Christoffel transformations. They map the upper half of the complex i-plane 
onto the interior of a polygon in such a way that the real axis is mapped onto 
the boundary. 

Referring to FiglJl the mapping is 

^ = 5(t - a)("A)-i(t - 6)W-)-i(t _ c)(-^/-)-i . (1) 

where S" is a constant of scale and a, 6, c, • • • are the point of the real axis corre- 
sponding to the polygon vertices. 

If the polygon is a rectangle, it has four vertices and all angles equal to 7r/2. In 
this case EqUJcan be written in integral form as 

z{t) = S [ dt'{t' -a)-i{t' -b)-i{t' -c)-i{t' -d)-i (2) 

J to 

This 4 parameter class of integrals can be reduced to the following one 
parameter class to be discussed in the next section 

z{t,k) = S' di' [(l-i'2)(l-fc2<'2)] 5 (3) 

2.3. Elliptic integrals and functions 

A large class of integrals involving rational functions of root square of poly- 
nomials up to the fourth degree can be reduced to three types of parametrized 
integrals, the elliptic integrals of first, second and third kinds. 
The elliptic integral of first kind in its Legendre normal form is the expression 
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Figure 2; Example of Schwartz-Christoffel transformation 
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in Eq[3]with S' = I. For a general t, the integral is said to be incomplete and 
is given the symbol F{t,k), where k is called the modulus. When t = 1, it is 
complete and is given the symbol K(k) ~ F{1, k). Derived from k is the com- 



plementary modulus, k' = \/l — k'^. Another definition is K'{k) — K{k'). 
In H the following important relation is derived 



F{±llk,k)^±K{k)-iK{k') (4) 
Other relations deduced fr'om the definition of F{t, k) are 



F{t,k) - F{^,k) = -F{l,k)+F{^^,k) t>i 

F{t,k) = -iK{k')+F{^^,k) (5) 

EqlSl reduces the calculation of F{t, k) for t> ^ io the evaluation of F{t, k) 
for t <\. In particular limt_>±oo F{t, k) = —iK{k'). 

Inverting EqI3]with 5' = 1, t is expressed as function of z and k, defining the 
function sine amplitude sn 



t = sn{z, k) (6) 

2.4- Capacitance calculation with Schwartz- Christoff el mapping 

The application of conformal mapping to capacitance calculation in geom- 
etry with 2D symmetry is based on the fact that capacitance is preserved by 
conformal transformations. For regions with strips and planes on the boundary, 
the mapping is a Schwartz-Christoffel transformation 

The basic configuration is shown in FigO where the upper half of the t-plane 
is mapped into the interior of the rectangle in the z-plane. 
The mapping is expressed by elliptic integrals and functions 



^ - 2K{kf^Wt/r^^ ^ ~ Wt + 2St 

where Wt and St are the strip and gap widths in the i-plane. 

This function maps the strip in the t-plane onto the upper side of the rectangle 

in the z-plane and the two lateral half-planes onto the lower side. In the z-plane 
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the rectangle width is Wz = 1 and the height is = 2K(k) ■ 
The reverse mapping is 

t='^sn{2Kik)z,k) (8) 

The capacitance between the strip and the planes in the upper half of the i-plane 
is equal to the parallel plane capacitance in the z-plane 

where e is the dielectric constant of the material. 
If the upper and lower half-planes of the i-plane are filled with media with 
different dielectric constants ei and £2, the field component normal to the non- 
conducting boundaries is assumed to be zero, so that the calculation is done 
separately for the two half-planes. 

This approach has been pionereed in [l| where several structures approximating 
a micro-strip detector section were studied. The major limitation in the analyt- 
ical calculation using conformal transformation consisted in approximating the 
lateral strips as a single lateral half-plane so that C„ could not be estimated 
separately. 

This estimation was possible using a mixed analytical- numerical method Qj-Q 
that, beyond conformal transformation, involves numerical integration and ma- 
trix inversion. 

In this paper this limitation is removed and the relevant capacitances are ex- 
pressed by analytical expressions. 

2.5. Auxiliary structure 

The structure in Figgis used in the following as an intermediate configura- 
tion and we need the capacitance of the strip to the opposite plane. 
This capacitance is calculated applying two transformations. The first one maps 
the interior of the rectangle in the z plane onto the upper half of the i-plane 



(10) 



where the strip width Wz and the left and right half-gaps Sl and SI are 
defined in FigUl The modulus ki is defined by 
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This transcendental equation can be approximately, but with high precision, 
solved employing a series expansion from j6||. The first order is 



K'(k) _ ^M2i±^) 0<fc<^ 

I l„(2i±^) V2- - 

that can be inverted to obtain k. defining y — as 



^, ^ / exp-^-2. y :^>i 



(12) 



Vexp^?^+2./ i4:(fc) 
/exp'^/j^-2.\^ K'{k) 



\exp''/y+2.J K(k) 
In the i-planc the strip width and the half-gaps are 



< 1 (13) 



With the definitions 



the mapping from the upper half of the t-plane to the interior of the rectangle 
in the w-plane is 

, , /■* dt' 
u{t) = / , (16) 



s[-w,- s: 



Following the steps outlined in [J|-[lO|-[ll|, the capacitance is 



K'{k,) Sj SI 

^-'K-'ii{k^ ~ y mT^y w;tsi ^^^^ 
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t— plane 



(1) 



(2) 



W =1 



(4^ 



z-plane 



Figure 3: Mapping strip with two planes into a rectangular domain. Top: original structure, 
upper half-plane; bottom: parallel plane capacitor 
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3. Strip to ground capacitance 



The strip to ground capacitance Cg is calculated applying the Schwartz- 
Christoffel mapping from the rectangular domain in the z-plane to the upper 
half of the t-plane and then the reverse mapping to the rectangular domain in 
the u-plane with a different modulus as in FigO 
The mapping from z plane to t plane is 




(18) 



where the modulus ki is obtained implicitly from 



(19) 



that can be solved using Eq fT^lfT51 

The mapping from t plane to u plane is the elliptic integral 




(20) 



where fc2 = kisn{K{ki) 
The capacitance is 



, ki) and ki come from Eq[Tn]or EqfT^ 




2K{k2) 

K'{k2) 



(21) 



In FiglH Cg normalized to the parallel plate capacitance Cj 
plotted versus the yj^^g with ^'f^^' as parameter. 



'pp ~ 
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Figure 4: Mapping sequence for calculating C in a rectangle with asymmetric gaps. Top: 
original structure; middle: intermediate structure, upper half-plane; bottom: parallel plane 
capacitor 
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Figured: Mapping sequence for calculating Cg . Top: original structure; middle: intermediate 
structure, upper half-plane; bottom: parallel plane capacitor 
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w/(w+s) 

strip to ground relative capacitance versus w/(w+s) 



Figure 6: Cg normalized to parallel plane capacitance Cpp versus and "'^^ {w is the 

strip width, s the interstrip gap and h the height). The curves with symbols are the results 
of the calculations in this paper. The curves without symbols are the results of the empirical 
model proposed in 12]. 
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4. Interstrip capacitances 

4-.1. Preliminary considerations 

The geometry of the detector relevant for the study of interstrip capacitances 
is displayed in FiglT] The upper half of the t plane is in air, the lower half in 
the detector material, silicon in the following examples. 

A simplified assumption is that the silicon section extends indefinitely, that is 
h oo. In this case each interstrip capacitance is obtained multiplying the 
half-plane capacitance by 1 + e^' (e^' is the relative silicon dielectric constant). 
For later use we assume that the point on each gap where the electric field 
generated by a voltage applied on the central strip is normal to the detector 
surface is in the gap center (interstrip points). These points must exist for 
continuity of the electric field and the field lines connecting them with the 
central strip divide the detector in separate regions. For each lateral strip, 
these regions are the volumes where a test positive charge would drift under the 
action of the electric field due to the central strip. 

In FigIT] the coordinates of the strip edges Xm , Xon and the gap centers Xcn of 
a small section of the detector are labeled. The interstrip capacitances in the 
upper half-plane are estimated with the help of FiglH Xcn are the interstrip 
points. 

The points marked on the picture are: 

x^n = niWt + St)-Wt/2 (22) 

Xon - niWt + St) + Wt/2 (23) 

Xcn = {n + l/2){Wt + St) (24) 

The reason for introducing these interstrip points is that it can be proved 
that in the rectangular domain the field lines running from these points to 
the opposite plane are straight lines normal to the planes. Those divide the 
rectangular domain in smaller rectangles with the geometry shown in Figd] in 
which the Poisson equation can be solved separately. 

4-2. Capacitances between neighbouring and not-neighbouring strips 

In FiglHl the upper half of the f-plane is mapped onto the interior of the 
rectangle in the u-plane by the transformation 

1 t . \ . Wt 
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The central strip in the t plane is mapped onto the upper side in the u plane 
with width Wu ~ 2. The gaps on the side of the central strip are mapped onto 
the lateral sides of the rectangle. 
EqHl gives 

hu = K{k[) (26) 

and Eql5] implies that all lateral strips are mapped onto strips on the lower 
side of the rectangle. 

The n — th lateral strip is mapped onto a strip with width 

Wun = F{Xon/{Wt/2),h)^F{x,„/{Wt/2),ki) 

ki Wt ki Wt 

The left and right gaps of the n — th lateral strip are mapped onto gaps with 
widths (defining n = 1 F{xcn-i/{Wt/2), fci)) = 0) 



SL = F{Xcn/{Wt/2),ki)^F{Xon/{Wt/2),ki) 

Wf + St 1 Wf + Sf 

= F((2n+l)^^,A;i)-F(-- + (2n+l)^^,A;i) (28) 



= F{x,n/{Wt/2),ki)-F{x,^^i/{Wt/2),ki) 

1 Wf + Sf Wf + Sf 

= F{- + {2n~l){^^),k,)-F{{2n+l)^^,k,) (29) 

The problem of calculating C„ in a half-plane is reduced to the calculation of 
the strip to backplane capacitance in a structure like FigUl following the steps 
described previously. The contributions in air and in silicon are added. 
The steps are the foUowung: a modulus fci is determined by 

Wun + SL + SL o Kik 



2-^^ (30) 
K'{k,) ^^^^ 



using EqdSl EqfT4l gives 
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Si 



n 



) 



(33) 



(32) 



(34) 



and therefore the n — th order interstrip capacitance is 



Cn = e 




Ql I Or 



(35) 



The interstrip capacitance Ci between neighboring strips is displayed in FiglHl 
versus W/{W + S). The higher order interstrip capacitances C„ for n E (2, 7) 
are displayed in Fig[Tni Ci depends strongly on W/{W + S) while the C„ de- 
pendence is much weaker. 

Fig ll II displays the total interstrip capacitance Cis^tot obtained adding all inter- 
strip capacitances up to the 7-th order on both sides and the total capacitance 
Ctot = Cis,tot + Cg. These quantities are easier to measure than C„ and to 
compare with experimental data. 



5. Comparison with other approaches and experimental results 

The results presented so far refer to geometrical capacitances in fully de- 
pleted detectors with strips on the junction side. Realistic detectors may differ 
from this model for various reasons. The most relevant is the presence of free 
charges below the interstrip oxide acting as additional conducting layers. This 
effect is difficult to predict and describe and may alter significantly the mea- 
surements. 

For detectors not satisfying the previous conditions, either our model is made 
more realistic, and more difficult to solve analytically, or a numerical approach 
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Figure 7: Section of a micro-strip detector with strip edges and interstrip gap centers marked. 
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t-plane 



W, % St W, s, Wt St W, S, W, W, 



'^ol ^^12 ^o2 ^c2 ^oi 

k,-wy(w.+2so F(t/(wy2),k,) 



u— plane 




s^(2K(k,)z/W„,kY)wy2 



Figure 8: Mapping of the half-plane in air. Top: original structure with few strips, the strips 
extends indefinitely on both sides. Bottom: rectangular domain with all lateral strips mapped 
on the lower side. 
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Figure 10: Plot of the higher order interstrip capacitances Cn versus w/{w + s) (w is the strip 
width, s the interstrip gap and h the height). 
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Figure 11: Plot of the total interstrip Cis^tot and total Ctot capacitances versus w/{w + s) for 
w + s = bOfim and h = 300/im. [w is the strip width, s the interstrip gap and h the height). 
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like that in 13] can complement our analytical approach by studying the effects 
of deviation from ideality. 

In the following, the field of applicability of the calculations developed in this 
paper is studied b y co mparing them with other models and the experimental re- 



sults presented in 
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14{ . These results have been selected because they refer to 



fully depleted detectors with strips on the junction side and scan a broad range 
of detector parameters. Furthermore data are presented in a format adequate 
for comparison with calculations. 

5.1. Cg comparison 

In Figl6]the values of Cg calculated with our model are compared with the 
empirical formula developed in [l^ obtained fitting the results of the numerical 
solution of the Poisson equation presented in [isl l . This formula is expressed as 

(36) 

—h-J\w+s> 

where /( ^^g ) is defined in [l^ . This model is supported by experimental 
results presented in 

In [l^ results for different material orientations (111) and (100) are presented. 
(Ill) has an higher trapped charge at the Si-02 interface and is more sensitive 
to irradiation. Nevertheless for nonirradiated devices, in the frequency range 
relevant for detector signal and for adequately high overdepletion voltage the 
capacitances are not dependent on the material orientation but only on the 
detector geometry and can be compared to our calculations. 
FigElshows an excellent agreement between our results and Eq[3ni Therefore, at 
least for Cg , detectors can be operated in such a condition that the geometrical 
capacitances are experimentally measurable and they are described very well by 
our calculations. 

5.2. Cn comparison 

A comparison of C„ is possible with the model developed in [ll| , where the 
Cn for n E (2, 7) were calculated for W ~ S. The results are compared in Tab[T] 
showing an excellent agreement 

Measuring separately C„ is quite difficult. The most straightforward compar- 
ison is with the total capacitance Ctot (including Cg) or with the total interstrip 
capacitance Cis^tot- 
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Table 1: Cn (fF /cm) calculated in and in this paper for W = S 





Ci 


C2 


C3 


C4 


C5 




C7 




487. 


96. 


40. 


22. 


13. 


9. 


6. 




478. 


95. 


39. 


21. 


13. 


9. 


7. 



5.5. Ciot comparison 

Experimental measurements of Ctot are presented in In order to 

compare the results, it is necessary to fix also {W + S)/h. Ctot and Cis^tot 
shown in Fig[TT] are calculated with W + S ^ 5Q^m and h = 300fim. The 
easiest comparison is with the formula, supported by measurements on fully 
depleted detectors with strips on the junction side, presented in [ij], 

w 

Ctot = 0.8 + 1.6 pF/cm (37) 

w + s 

A linear fit to Ctot in Fig[TT] (not shown) gives 

w 

Ctot = 0.73 + 1.60 pF/cm (38) 

w + s 

showing an excellent agreement. 
6. Conclusions 

The capacitances in micro-strip detectors are calculated analytically with 
the help of Schwartz-Christoffcl conformal transformations. The results are ex- 
pressed with elliptic functions. 

The strip to ground, the total and the interstrip capacitances are calculated as 
function of the geometry of the detector. 

The results are compared with other models and with empirical formulae sup- 
ported by experimental measurements showing an excellent agreement for non- 
irradiated and fully depleted detectors. 

Appendix 

The numerical evaluation of the complete and incomplete elliptic integrals 
and of the elliptic functions is based on FORTRAN 77 routines from the CERN 
library package available in 11611 . These routines are based on numerical algo- 



rithms developed in 



lable in hj. T 

0Q,liIQ- 
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